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AHAI<YS1‘G AlID APPWCATIOH OP rarilMUH VARIANCE 
DinCRF.TR BYSTEM IDPNTIPT CATION^ 

S. ICotob^^** and' K. JCauTmn*** 

Abstrac t 

■ ptiinal design oT an on-llno state and parameter estiiration algorithm 
results in a complex structuj-G not readily feasible for adaptive control pur- 
poses and/or for Implementation in a typical process control or flight computer. 
Suboptimal designs are often used in practice and despite the ease of imple- 
menting these procedures, problemB such a.s di 'erpence and/or inacciu’acie.s are 
not unconunon. An on-line minimuin variance parameter identifier has therefore 
been developed which embodies both accuracy and computati nal efficiency. The 
new forinu!i ,tion results in a linear estimation problem with both additive and 

t 

multiplicative noise, 'ilio resulting filter is shown to utili^.e both the co- ^ 
variance of the parameter vector itself and the covariance of the errc>' in 
identification. 

A bias reduction scheme can be used if desired to yield asymptotically 
unbiased estimates. It is proven that the identification filter is m.ean sa.uare 
coirvergent and mean square consistent. Ibe If/ parameter identification scheme 
is then used to construct a stable state and parameter estimation algorithm. 
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1. INTRODUCTION 

Simple mechanical linkages are often unable to cope with many control problems 

associated with high performance aircraft. This has let to the present interest 

fl 2l 

in digital fly-by-vire ’ control systems. Digital implementation is extremely 
advantageous because of{ 

» The significant weight and vol\une savings. 

• The availability of integrated circuits. 

'' The ability to design complex controllers vdiloh were previously 
impossible to implement on board an aircraft. 

• The high reliability of digital logic. 

Furthermore, the need for an adaptive control system has been established for:^^’^*^ 

• providing uniform stability and handling qualities over the 
complete flight envelope despite drastic changes in the open 
loop characteristics of the aircraft. 

• providing, acceptable flying qualities over a wide range of 
external disturbances due to atmospheric tui'bulence and 
outer loop command signals. 

In designing an adaptive control system, it is necessary to determine whether 
to implement an explicit system in which on-line pai'ameter identification is per- 
formed or an implicit system which does not require explicit parameter identifica- 
tion. Recent studies have indicated preference for explicit designs whenever the 
process to be controlled has non-minimum phase characteristics and/or high gain 

[51 

and large bandwidth limitations , 

Tlie development of an implementable digital adaptive control system requires 
the use of an identification schume that is capable of supplying parameter esti- 
mates at an accuracy and rate specified by the controller characteristics. Because 
a digital adaptive controller uses elements of the discretized matrices, identifi- 
cation of these elements and not the continuous physical system parameter should 
be considered. Furthermore, identification of the parameters of a continuous 


system (e.g., stability derivatives) from discrete data results in a problem with 
many severe nonlinearities, 

Frequently* the chotoe of an identification method depends not only upon 
the type of model used* but also on the computinc facilities available. Duo to 
inherent limitations in the storn(;o and computation time of a typical process 
control computer the identification scheme should be recursive in nature In 
order to avoid data accumulation. Therefore, of interest were the extended 
Kalman filter, the weighted least squares algorithm, stochastic approximation, 
and a decoupling process in which linear state estimation and parameter identifi- 
cation were considered separately and alternate] y. 

Despl ' he case of implementing these procedures, problems such as divergence 
and/or inacciu’acies are not uncommon. In particular the following items were 
noted: 

♦ Tlie pei'formanco of the extended Kalman filter is inadequate 
for adaijtive control purposes because of inaccuracies wdiioh 
result from linearization^ . 

• The Weighted least square procedure is very useful for high 
signal to noise ratios. Highly biased estimates are not un- 
common at low signal to noise ratios and in the estiiriates of 
insensitive parameters^ . 

• Computer sijnulatlons conducted using the stochastic approxi- 
mation (RA) algorithm for the identification of linear lateral 
motion of an aircraft showed that the convergence of the SA 

[ 12 ] 

procedure is very slow especially for adaptive control purposes. 

* A decoupling process in which linear state and parameter esti- 
mation are performed separately and alternately ( state estimates 
used in the parameter filter end parameter estimates used in the 
state filter) is highly Inaccurate and of limited use in digital 
adaptive control. 
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Linear ayfatem iclontificatlon using the input and noisy measurements 
of the output can be generally oast as n state estimation probleiti vitli both 
additive and multiplication noise (A^4N). 'Ilieso terms will in fact be functions 
of the same noise sequence. The continuous optimal nonlinear filter ns derived 
by Kushner^ for (AI»!N) is infinite dimensional and its physical realisiation is 
impossible. Approximate linear filters were subsequently derived for 
under che assumption that the additive and multiplication disturbance' terms are 
functions of two independent i*andom processes; hence those results are not 
immediately applicable to system identification. Thus a new on-line mlnimuin 
variance filter for the identification of systems with additive and multiplica- 
tive noise has been developed which embodies both aocurncy and computational 
efficiency. T5ie resulting filter is shown to utiliae both the covariance of 
the parameter vector itself and the covariance of the error in identification. 

A bins rednction scheme can be used if desired, to yield asymptotically imbiased 
estimates. 

As common in deriving anj'' estima.tion scheme, proof of the convergence 
of the identification filter is an integral part of the validation of the result 
In this respect, the proposed identification scheme is shown to be convergent In 
the mean square sense. Tlie proof consists of deriving a suitable upper-boxind 
for the mean square error (MSE) and showing that the MSE converges to zero as 
time tends to infinity. The mean square convergence of ‘the filter implies con- 
vergence 'with probability Which, in turn, would imply that the estimates are con 
sistent. Using the proposed pai’ameter identification filter and the related 
convergence proofs, a state-parameter estimation scheme is consti’ucted and 
proven to be stable in the sense of boundness. Tiie resulting state-parameter 
scheme is shown to be computationally feasible and amenable for on-line system 
identification and adaptive control applications. 
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To Illustrate the reliability of the idenflcatlon schemes and the problems 


encountered, experimental results for simulated linearised lateral aircraft 


motion in a digital closed loop mode, are included. A comparison of the extend- 


ed Kalman filter and tlie minimum variance filter in the adaptive mode are presented. 


2. PROBLEM DEPiniTION 


The problem of detei’inlning on-line values of certain parameters apiJearlng in 


the discrete. eq.uationa of a lineal* constant coefficient process, VliicJi are best 


with regard to use in adaptive control logic, given the input and noisy measure- 


ments of the output, vas considered using an on-line minimum variance filter. 


'liie corrorponding equations are: 


x(k+l)=A x(k) + B u(k) 


= x(k) + n(k) 


vhere 


V il 

x(k) ^ plant state at the k stunple instant (n x l) 


A = state transition irntrix for the discrete system (n x n) 


u(k) “ Cf trol vector (m x 1) 


B “ Control distribution matrix (n x m) 

y(k) « measurement vector x 1) at the k^^^ instant 

T](k) = measurement noise at the k^^' instant lovariance 

matrix R(1,5) = O 5.,, 


.where “ 1» i “ ’^ij “ ^ otherwise. 


3. PARAMETER MODELING 


To e.3timate any unknovm vector of parameters q appearing in the .state trans- 


ition matrix A and in the control distribution matrix B, it is necessary to model 


the dj^namics and observations of the system parameters. Furthermore, since not 


all parameters appearing in the matrices A and B are to be identified, the differ- 


entiation between the set of parameters that are to be identified and the set of 


REPRODUCIBILITY OF THE 


APTAlTXTAt PAnm T.Q DA AD 


parametors not to Iiq identified is generally rocoiraaendcd. In particular, a 
convenient representation of the systen is; 

x(k+l) = C(k) q(k) + D(k) a (3) 


vhere 


C and D are selection jr.atrice3 containing values of the system 
otatc and control at the k instant. A aero entry for a part- 
icular C^j(or vould indicate that no coupling exists "between 
and (or between and s^), 

q is a vector of unknown pai'ar.eters appearing in the A and B matrices 
S is a vector of known parameters appearing in the A and B matrices, 
llio model for thq constant deterministic or stochastic parameter vector q is then 
given by: 


Systems Ih^nairdcs 

q(k+l) = q(k) E(q(0))=qQ (q(O)-qy ) (q(O)-qQ)'^] (»f) 

Observation 

x(k+l) - C(k) q(ii) + D(k) S (5) 

liie identification problem as defined in (h) and (5) appears to be a conventional 
linear state estimation problem. However, because x is not known exactly, C and 
D are also unknown, and the rules and usage of the conventional Kalman filter 
can not be applied. Substituting equation (s) into (5) gives: 

y(k+l)-n(k+l)‘=C(y(k}-ri(k) , u(k))q(k) + D(y(k)-ri(k) , u(k)} S (6a) 

where 

the notation C(y(k)-ri(k) , u(k)) and D(y(k)-n(k) ,u(k) ) stresses the fact 
that the selection matrices C and D are functions of the state and 
control values. Noting that: 

C(y(k)-n(k), u(k))=C(y(k), u(k))-c(n(k), u(k)) 
and similarly for D, equation 6a can be rewritten as: 
y(k+l)-D(y(k), u(k))S=c(y(k), u(k))g(k) 

-G(n(k), u(k))q(k)-D(r|(k), u(k))S+n(k+l) (6b) 
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Do fining 


z(k) = y(k+l) “ D(y(k), u(k)) S 

to be a pseuAo-mousuroinont vector for the linear oystfin given in {10 und rearrang 
Ing equation (7) givea the parameter ’''Xj 

q(k+l} « q(k) (10 

z(k)sC(k)q(k)-C^(k;q{k)+n{k+l)~D^(k)S (8) 

where 

C(k) = C(y(k),uO')) 

C^(k)=C(n(k)Vu{k)) 

r.^(h)=ii(n(k),u(k)) 

Equations (ll) and (8) denote a linear time invariant system with a transition 
TiiatriJ: I and observation matrix C(k). The observation is corrupted by the 
multiplicative noise term C^(k) and the additive noise tei’ins n(k)~D^(k)S with 
covariance 11 cp,+li{D S D '^) 

GO ^ n n 

h . MIIlimM VARIMCE ESTIliA.^^^ 

1<*1 Develrpmcnt 

Becaiise of divergence and/or inaccuracies cownon to most existing identi- 
fication schemes, it was desirable to develop an alternate scheme that could 
hopeliillj, deal '^ith these problems. The proposed filter is based on a minimum 
variance performance index for the state estimation of a linear system with 
additive and multiplicative noise. . 

Tl;e optimum minimum variance filter for a continuous system is in fact non- 
f7l 

linear , and its exact implementation is virtually impossible. A linear 
optimal filter was therefore of interest. Thus defining the identification 
algorithm to be 

q(k)-q(k-Ji)+K(k)[z(k)-C(k-l)q(k-;i)] (9) 
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K is to be detomlnod co as to minimize: 

J « i Ey(q^(H) - (10) 

Vhich is the trace of the covariance matrix 

P(k)=E J(q(k)-q(k))(a{k)»a(k))'^’} ( 11 ) 

M 

Where E„(a(k) )=S[ajy(0) , . . ,y(k) ] 

'ihe parameter A (\diich vrill bo discussed in the next section) defines the fre~ 
quency of identification, Define the parameter error as q(k)*-c(k)-q(k) , with 
Initial conditions at '1-0 

E(q(0))n0 ; E{q(0) q(o)^'} =P^- 

Usln/' (l|), (6) and (9) the error q propaf'iates as; 


q(k)=q(k-5-/H-K(k)[-.C(k-l)q(k-il)~C^k-l)q(k)-D^(k)r>+n(k)] _(l2) 

By post inultiplyinf/ 12 by its transpose , and taking the conditional expectation 
over the entire measurement vector history (y(o) .y(k) ) and noting that; 

E[C^(k-l)q(k)q^(k-A) |y^, . . .y(k) ]=0 

when Jl > 1, the diffei’ence equation for the concUt’onal variance P(k) becomes! 

p ( k )=P( k-5- ) -K( k ) C(k‘-1 ) P( k- 1 ) -P^’( k-il ) C^( k-1 ) K^(k ) 

+K(k)(S(k-l)P(k~A)C^’(k-l)'iw(k*-l)+Bg^)K''^(k) (l3a) 

vrhei’e rp 

W. =E{0 (k-l)q(k-;,)q'®(k-A)C;!;(k-l)} (13b) 

X.—J. n >1 

stationary conditions for the minimization of the trace of P(k) are obtained by 
setting all derivatives of (l2) with respect to the elements of K(k) equal to 
zero. This yields: 

K(k)=P(k-J!t)C^’(k-l)[C(k-l)P(k-A)c''^‘(k-l)-Hw(k-l)+E 1“^ (lij) 

k.2 Observations 

(l) Tlie gain of the resulting filter is a fxmiction of the error 
covariaiaee P and the weighted noise covariance w(l3b); where 
the weightinc; matrix for d) is the covariance of the identified 
parameter q. 


( 2 ) Ihe derivation of the proposed minimum variance filter is made 
possible by not identifying every sample; i.e., ^ > 1. For 
!, = 1, the expected value of many cross terms involving the 
parameter q» the error q and the noise selection matrix C„ will 
not vanish; this can be illustrated by noting that: 

EEPRODUCIBILITY OF THE 
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E{c^(k-l)n(k)a^'(k-l)|yQ,...y(k)} 0 

since Cp{k-l) and 5 (k-l) are correlated. Obviously rocuroive 
updating of non-vanishinp crosn terms is possible, and the 
minlmlaation vo\fLd not be feasible for on-line implementation. 

£!ittce no estimation is performed botver . the time instants 
k-jt and k-1, the notation used Instead of 

Hence the resulting identification algorith' can be given by 
the folloving equations: 

q( k )=q( k-i )i ''■( k ) [ s ( k ) -C{ k-1 ) q( k-fi.) 1 


K(k ) =P( k-1 )c''’(k-l ) [C(k-1 )P ( k-1 ) C''^(k-1 ) k-1 

P(k)=r(k-l)-K(k)C(k-l)?(k-l) 

vhere equation (l?) results from substitutinp. equation (lit) in (l3) 



!). nun REDUcnoH 


Althounh the linear jnlnimujn varlrmce filter as deacrihed hy oquationo ( 9 )> (l^t) I 
and (15) vaa obcorved to be relatively accurate with respect to other linear schemes, 

\ 

a substantial bias did appear in the paramete*' sstimatea, especially in the estimates | 
of insensitive parameters'^^. An investigation was therefore condtic tod to dotor- 
rnine the causes and the means to reduce or eliminate the blaG. 

Ey combining eqviations (15) and (it), the gain can be revidtten asj 
K{k)=P(h)C”(k-l)R“^(k-l) 

By substituting (16) in (9) and talcing the expectation and the limit ns K **■ <“, 
it is found that; 

E(q(lc)}=[lHC^(k-l)R"^- 6(k-l))+E(c'^'(k-l)R“^ 

C^(k-a)}]"^- • [R{c'^’(k»l)R^’^&(k-l)}] • q(k) (IT) 

Obviously equation ( 17 ) reveals the bias in the etJtiinetes of q, q. Assuming that 

the term E{C‘^(k“l)R^^C(k-l) } is a generalised measure of the signal power, and 

R{C^(k“l)R^ Cjj(k~l)}ls a generalized measure of the noise, equation ( 17 ) can be 

Written as: / . 

E{q(k)} = [S + N]~-^ rn] • q(k) (I8) 

Where 

S = E{C-^(k-l)R“-^ C(k-l)} 

N = E{C^(k-l)R“^ ^ 

By examining equation (l 7 )and (iQ), it becomes obvious that the troublesome term 
is the noise power N. Hence a correction term must be added so as to compensate 
for the bias. From consideration of equations (16) and (17), it is clear that the 
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corroctlon tsrm nmat incorporate the 
AdOlnc the correction term to (iH), 
q(k)«(l+P(k) G(U)}q(h-Jl) + 
where G{k) is to he found such that 

llirt E{q(K)} = 


covariance tovjti Pj^ and the latoat oBtimato 
the basic algorithm bocomesj 
K(h) [z(k)4{k-D 


a 


K'>“ 


( 19 ) 

(PO) 


Tidting the expectation for (l?) and using (20), yields! 

G(k) B E{C^(k-l)H“^Cj^(k-l)}«If 
Hence, the modified minimiun vai’iance filter is given by; 

q{k)=(l+P(k)JJ{cJ|(k--l)R“^- 'C^(k-l)})$(k-.A)+K(k)[z(k)-e(k-l)q(k-il)] (2l) 

where P(k) and K(k) are given by the recursive equations (l4) and (15). It should 
bo pointed out that in recursive on-line parameter identification schemoa, onlv 
asymptotic unbiasness is possible 

6. PIb:’r.rSf ABILITY 


Essential to any estimation scheme is the validity of the resulting ostimatos. 


In this respect, it is desired to prove that the proposed identification algorithm 


convf i,es to the actual system parairoters. The convergence of the filter is of 
particular iinportance since the resiAlting estimates are to be used in the consti’uct- 
don of OJ) adaptive controller. 


The proof of convergence of any Gstimation scheme usually involves the following; 

1. Proving that the natural modes of the estimation scheme 
are stable. 

2. Establishing, in a statistical sense that the errors 
resulting from the measurement and/or plant noise 
remain bmmded and/or converge to zero in the limit. 

Before proceeding in establishing the convergence of the proposed identification 

scheme, the following assunptions needed for the proof are stated: 

Al. is a vector sequence whose entries are zero mean independent 

variables . All entries of the measurement noise vector (fj,) are 
mutually independent. Second and fourth moments are uniformly 

bounded. ■ 
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A2. TJio flotei'mlniotic conbrol voctoi’ u(k) ia auumod to be bounded. 
Similurly the output vector y(k) and all trarml’ornationo on y(k) 
are uuniwod to have bounded momcnto. Clearly, thla aoBujnption 
ia roaliat'f ”• cinco the proposed alKorithm is to be used In idonti- 
lyinp phyca.v,ttlly I’calli'.ablc oyatoms, 


A3. I'ho linear syotoin is completely contro.t lublo uijd coiiipletoly 
observable. These conditions ere necessary since the only 
part of the oyatom vhich is identifiable from input~outpuk 
observations la the coinplotcly conti’ollablo and co»np3c'toly 
observablo one. 


All. 


The paraireter.net to be identified is assumed to be completely 
observable ^ ^ in the senno that the Information matrix 

'7J(k,D = I c'^(i)R“^ C(i) 
iol “ 


is positive definite. Furthormoro, it is assujned that tlie 
system identification model is uniformly completely obsorv- 
able in the sense that 

0<oI<c'^'(i)R7^ C(i) < Vi, a.s. 

u) ' 

vdiere a and 3 are constants. 


A5, Qhc product of the m.atilx P and the sinnal pov/er B is positive 
definite. 


6.1 ParnmetCr filter converreiice 

Tlieorcn 1 . I'oan Square Converpongo . Under the assumptions A1 to A5, the 
linear estimator of q(k) nlven in equations ll*, 15 tu^d SI conver/^es in the 
mean square sense to the unknown partuneter vector q of the linear system 
in (1-2). 

Proof ; 'Ihe estiination error q can be piven by tlio follovdnc equation; 
q( k ) = [ I+P k-1 ) K“^C ( k-1 ) ] q( k-p ) +P [ C-'^( k-1 ) n”^G( k-1 ) q~Nq 

-e'^(k-l)R“^C(k-l)q-c''^V.-l)P~^(n~D^S)] (22) 

Preinultlplyine (22) by q(k), takinp the exception and by repeated use of 
the Cauchy-Shwertz and triangle inequalities, an upper bound for the moon 
square of the identification error can be given by: 

E( I laOOl I =’•) < {1+1 |P| I =’J£5-2X E( I |q(k-A) I j 2) +1 |P| I ^^'6 (23) 

where is the minimum eigenvalue of the product PS, arid ^5 and ;C6 are 

nonliiisar functions of the signal and noise power, 


U 


The upperhound in (23) ia ti /'onorullijution of tho bovuid Riven by Mondol^^'®^ 

for ulnRlQ output Idcntiflsation. A dokuiled derivation of (23) io provided 

fill 

in roferencQ 12. By upplyinR Venter 'n 'i’Jjeorejn ^ and usiiiR tho fact that in 

1 

the limit, the rnaximujn elRcnvaluo of P behavoo ao — , the mean cauaro error 
B[ I Ifi(h) I j is Dhovn to conyerRG to zero in tho limit, i.e., 

Aim 1![ I !<i(k) 1 1 “ 0 . 

k-»w 

Combining the fact that the proposed identification achemo is asymptotically 
unbiased and mean SQuare convorRont, it is concluded that the filter is mean'- 

f « I h 

square conaiotent^“ . 

6.2 ConverReneo of f»tate Kstiratlon 

'ibe maximum likelihood, minir.iun variance and leant squares estimate of the 

state vector x(k) Riven the reasurement vector y(0},...y(k) is riven by the Kalman- 

ritl 

Bucy filter . 'Xhe Kaljaan filter was ohow to converpe in the m.ean square conKe 
and vith probability 1 if tho plant mode] and Gaussian noise statistics are 
exactly knov.'n. T- .,*3 vhere tho plant model Is not exactly knovm, an approx- 
imate Kalman filter can be constructed using identified pnrameter.s. Tlio stability 
of the approximate Kalman filter is diseuseod In the sequel. 

Thonrem 2: Given tlie approximate Kalman filter 

x(kyk)=A x(k-l/k-l)+B u(k-l)+K (k) y(k)-fA x(k-l/k'-l)+B u(k-l)] 

® - 

Kg(k)=P^(k)B”^ (2l() 

vhere state ertlmate using identified parameters 

A , B identified system and input matrix 

A' " . A 

and covariance matrix using identified parameters, 

if the linear system (l,2) is stable and if A and B are consistent estimates of 
A and B respectively, then 

n^ n ^ |x(kjk) - x(kjk)|=0 vith probability one 

nmOBUOIBILlTY OF THE 

original rage is poor 
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P (k) P° as k •' " vit'i) probability one 

D O 

K (k) ->• K® as k *+• w vith probability one 
vhere ° ^ 

x(k|k) ^ optimal state estimate 

K° , P° optimal /^ain and covariance matrix 
The proof of theorem 2 is provided in reference (lil). 

It vras shown in section 6.1 that the estimator q, in (21) is mean square 
convei'Korit and moan square consistent. By invoking theorem 2 ond the stability 
condition on the plant in 1, it can be observed that the proposed identifier 

t 

and the respective approximate Kalman filter constitute a stable (in the sense 
of boundedx)OS3 ) stutc-iiaruineter identification scheme. This structure can serve 
as an alternative to the lineariaed Kalman filter vith the advantages of stability 
and ease of implementation. Riinv lations for the state parameter filter’ are pre- 
sented in section 8. 

Y. IDENTIFICATION OF TIIIE VAHYUfo PARAiKTFRS 

An linportcmt property of p. parameter idenfication procedure is its effect! vu- 
ness in identiiying and tracking time varying parameters. The extended Kalman 
filter was of limited use due to problems related to lineariaation and the up- 
dating of the priming trajectory Furthermore, because of problems related 
to the convergence speed and tipdating of the correction gain, stochastic 
approximation was found to be unsatisfactory for identification of time vary- 
ing parameters. 

Modeling time varying parameters as a first order random walk, the minimum 
variance filter esn be modified or rederivod so as to track the variation in the 
system parameter q( It ) . The new parameter model is then given by the folioving 
equations : 
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Paratneter Itoclol 


q(k) « q(k“A) + v(k-^.) 

where v Is a zero meun xmcorrelated atutionary GauBsion noioe seciuencc 
with covariance ?i(k) » E(v(k) v(k)'*^). 

Observation 

z(k)=c(k~l)q(k)-C^(k)q(k)-Dj^(k)s(k) 

Minlyu.m Variance Filter 

q(k)=5q(k-P.)+K{k)tz{k)-.C(k-l)q(k-A)] (25) 

M(k) - P(k-l) + Q (26) 

K(k)=M(k)6'-^’(k~l)[c(k-i)ii(k}c'-^(k-l)+w(k-l)+R^J“^ (27) 

P(k)=M(k) - IC(k) C(k~l) M(k) (28) 

where (bCk-l) is the veif'hted noiae covariance. 

w(k-l)=E{C^(k-l) q(k-A) o.'-^’(k-A) cj*j'(k-i)} 

Since the parameter vector q is modeled as n first order random imlk, its 
covai’iance has to be updated recursively so as to compute w(k~l). Assuming that 
the initial parameter covariance T is given by; 

'f(O) = E{q(0) a7(0)}; 

the parameter covariance T(k) and the weighted noise covariance can be given by: 


T(k) = T(k~ji) +Q(k-A) 

(29) 

w(' .l)=E{C|^(k~l)T(k~.<l.)Cj^(k-Jl)} 

(30) 


Equations (25) to (29) summarize the minimum variance filter for identification of 
varying parameters. The resvU.ting filter is relatively straple for use in a typical 
pri-cass control computer. 

8. APPLICATIONS AND RESULTS 

The performanoG of the minimum variance filter was evaluated experimentally 
using a linearized model of the lateral motion of a typical fighter nircraft. 


Tliis evaluation consslsted of: 

(l) An inveBtigatlon of the convergence properties 
and parameter tracking of the minimum variance 
filter. 


(2) A study of the overall performance of the closed 
loop system in the adaptive mode. 

(3) A comparative EUialysis of the performance of the 

minimum variance filter: the extended Kalman 

filter and the weighted least squares algorithm. 

8.1 Adaptive Control Loop Design 

Development of an adaptive control system requires that consideration be 
given to designing a control .■'Igorithm that performs well and at the some time 
is easily adjustable on-line in response to parameter changes. 

Given a plant to be controlled and a model having the desired closed loop 
plant response, a model following adaiitive controller can be designed so that 
the closed loop system will behave in the some manner as the model, i.e., for a 
given input the closed loop system should respond ns the model would, were it 


subject to the same input. 

The particular controller used in the experiments was the one designed by 
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Alag and Kauiman- so as to yield a bounded error. In particu-lar for a plant 
of the form 

X (k+1) = A X (k) + B u (k) 

P P P P P 

and a model of the form 


x„(k+l) *= A„x_(k) + B u (k) , 
m mm mm 

the corresponding control signal u^ becomes: 


where 


u c u + u„ 
p 12 


U ={B^' B yh'^iA -A )x +(B^B )"^B^B u 
•L p p; p m p m p p' p m ti 

u„ = -K(x -X ) 

2 m p 


and K is such that (A ~E K) is stable. 

P p: ; ^ 
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Note that vhoreaD the feedforward cains that define eon be computed by 
formula evaluation, Given eatiinatcB of and Riccati typo iterative pro- 
cedure ic needed to find 
8 . 2 Computer Simulations 

In order to evaluate the overall performance, data for a typical aircraft 

was provided by NASA Langley Research Center for lineariiiod lateral airci*aft 

motions about trim^^^. !The physical aipnificance of the state and control 

vector components for aircraft lateral control is as follows: 

p \ roll rate 6 

u = a 

r yavr rate 16 

\ r 

8 sideslip anele 
4) I roll angle 

For simulation puj'poses , the aircraft was assumed to be flying in a fixed flight 
condition (FC2, Much 0.9, 3000 w) . A sensitivity study defined, for identification, 
a set of 12 parameters which make up the first and third rows of A and }3 matrices. 

To test out the minimum variance filter and adaptive controller in a realis- 
tic environment, the system shown in Figure 1 was simulated on an IBM 360/6? com- 
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pater. Parameter estimates were obtained using noisy measurements of the states 
Using Z=2, parameter ostlmtes vrere obtained alternatively every O.)^ sec. (sampling 
period =0,2 sec.) and were used every 1 see. in the gain adaptation procedure. 

The x’esulting parameters and gains vrere then used to estimate the states and con- 
trols every 0.2 see. The square wave aileron pilot input u^^ of 5° at the fre- 
quency of 0.4 Hz was used in all the experiments. 

The convergence properties, adaptive controller results and comparisons with 
different identification procedures to be presented in this paper were all con- 
ducted for the fixed flight condition (FC2) where all parameter estimates were 
initialized at 50 ^ of their actual values. 
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aileron deflection 
rudder deflection 



Fiifares 2 .a,b,c,d and e illustrate the behavior using the minimum variance 
filter, 'fhe asymptotic unblasuesa is evident in fi(f;ures 2 .a,c and e vhich shov 
that the estimates have converged vithin 20 sec. (50 measui’ements). As can be 
seen in figures 2 . b and d the estimates of parameters ^33 v<sre highly 

biased. The parameter ^‘eatured in figure 2 .b is very InsenBitive especi lly 
vith aileron ejccitation (the parameter* a^^ couples the sideslip angle to the roll 
rate ) . 

Figure 3 depicts the roll rate behavior in the adaptive mode. It can be 
seen that model folloving perfox'mance is highly correlated vith convergence of 
the parameter estimates. Reasonable model folloving vas achieved after 15 sec. 
when most parameters had converged to the actual values. 


To compare the perfonnance of the minimum variance filter, veightod least 
squares and the conventional extended Kiiiman filter (EKF), figures 2, t and 5 
can be used. The advantages of using the ItVF over the FICP is evident from 
figures 2 and 5. In v 7 oighted least fiquni‘c.s, onlj' 't pui'ameter estimates vere 
more accurate than those i*esulting from the IWF, (e.g. , as chovn in figures 2 .d 
and h. d); but the overall parameter estimation vms substantially better in MVF; 
where 7 parameter estimates vere more accui'ate than the estimates from the EKF 
and V/LS. 

It vas observed that estimation of insensitive parameters using M’/F vas 
relatively accurate compared to the other tested schemes. Since sensitivity 
is generally related to restrictions on types, pover and frequency of input 
signals as veil as the physical structure of the system; this method cou].d 
be very valuable in modeling and identification of complex systems such as a 
Fluid Catalytic Cracker at lov test signal povers and frequency. 
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9. DIBCURBIOnS ATfB CONCLUSIONS 

A linear minijnujii variunco pnvtimotier identifj.or vao derived and wan ahovn 
exporiinontally to eonvsrce to the actual pavainetcra. A bias veCuction Bcliemo 
and modiriontions i’or time varylnp paruiiioters wore preoented. The new filter 
pi’ovod aupcrior over oxintiiifi linear and linearlned purainetor filtorn ami 
conerally moi’e flc^dblo and effective in tlio oatiimition of inoonsitive paratnetero 
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